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ABSTRACT 

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present 
a new effort, with IRAM 30-m observations of ^^CO(l-O) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 
12-klog(0/H)~7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Z©); to 
our knowledge this is the largest number of ^^CO(l-O) detections ever reported for galaxies with 12-klog(0/H)^8 (0.2 Z©) outside the 
Local Group. We calculate stellar masses, Mgtar, and star-formation rates (SFRs), and analyze our results in conjunction with galaxy 
samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, T^ep, with Mstar and 
specific SFR (sSFR), we find a variation in r^ep of a factor of 200 or more (from ^50Myr to ~10 Gyr) over a spread of ~ 10^ in sSFR 
and Mstar- We exploit the variation of Tdep to constrain the CO-to-H 2 mass conversion factor aco at low metallicity, and assuming 
a power-law variation find aco^Z/Z©"^, similar to results based on dust continuum measurements compared with gas mass. By 
including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are 
metal poor, of low mass, and with high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases 
show that the models are generally quite successful, but at metallicities Z/Z©^0.2, there are some discrepancies. 
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1. Introduction 


Dense clumps in cool molecular clouds are thought to be the 
“cradles” of star formation, but finding molecules at low metal¬ 
licities has been exceedingly difficult. In dwarf galaxies with 
nebular oxygen abundances ^ 0.3 7© (12-hlog(0/H)^8.2Q de¬ 
spite vigo rous on-going star formation, there seems to be very 
litt le CO (|Sage et al.||1992[ [Taylor et al.||1998[ Gondhalekar et 
al. 1998t [Barone et al.||2QQ0[ Leroy et al.||20Q5t [Buyle et aL 


|Leroy et al.||2007[ |Schmba et al.||2()12[ [Cormier et al 


2006 


2014). Intriguingly, this is the same metallicity “threshold” be¬ 


low which little or no polycyclic aromatic h ydrocarbon (PAH) 
emission is detected in metal-poor galaxies ([Engelbracht et al. 


[20051 IWuet al.|2006[ [Madden et al.|20061|Hunt et al.|2010D 

The reason for this deficit in CO at low metallicities is not 
yet completely understood, although there is strong evidence for 


Send offprint requests to: L. K. Hunt 
^ Based on observations carried out with the IRAM 30m; 
supported by the INSU/CNRS (France), MPG (Germany), 
(Spain). 


IRAM is 
and IGN 


^ Throughout this work, we will rely on the |Asplund et al.| ( |2009 l 
solar abundance calibration of 12-klog(0/H) = 8.69. 


“CO-dark” molecular gas in the interstellar medium (ISM) (e.g.. 


Wolfire et al.|[2010t [Leroy et al.|[2011t Planck Collaboration et 


al.||2011[ [Paradis et al.[|2012| [Lebouteiller et al.||2012t [Pineda 
et al.[|2014 |. Such a component would result from a diffuse or 
metal-poor ISM in which CO suffers from enhanced photodis- 
sociation relative to H 2 (e.g., [Maloney & Black||1988[ [Bolatto 


et al.|| 19991 [Glover ^ Mac Low|[2011[ [Narayanan et al.|[2012|) 

H 2 can effectively avoid the photodissociation effects of UV ra¬ 
diation through self-shielding, even at relatively low gas surface 
densities, Xjsm- On the other hand, CO needs higher Xjsm to self¬ 
shield, thus implying that the CO-to-H 2 rati o should vary even 
within an individual molecular cloud (e.g., Rollig et al. 2Q06[ 
Leroy et al. [201 3 aj [Hughes et al.[2013 1. Under conditions of low 
^iSM and low dust extinction Ay, most of the carbon would be in 
other forms such as C i or Cii, rather than in CO. 


Nevertheless, even under extreme conditions in a metal-poor 
ISM, it is likely that stars form in cool, dense molecular clouds 
(e.g., |Krumholz et al.|2Qll| ). However, H 2 , the dominant molec¬ 
ular component of giant molecular clouds, is not directly observ¬ 
able, but rather inferred from line emission of rotational transi¬ 
tions of ^^CO(l-O), unfortunately deficient at low metallicity. 
Constraining the properties of metal-poor molecular clouds ob- 
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servationally is difficult because the factor, aco^ relating CO lu¬ 
minosity to H 2 mass Mh 2 can vary dramatically, because 
of the presence of CO-dark gas. aco can be more than 50 times 
higher than the Galactic value in metal-poor environments and 
in regions of low Ay and low l^isu (|Israel||19'97} |Leroy et al. 


2009bl[MTn|Pineda et al |20T0l [Bolatto et al.|20H t|Schruba et 


rich in atomic gas (see Sect. [^. This section describes particular 
features of the individual sources (in alphabetical order by name 
as in Table [T]). 


★ CGCG 007-025 
This galaxy is the 


very clumpy morphology as seen in Fig.[^ The radiation field is 
sufficiently hard to excite [Fev], a line that requires ionizing ra- 


diation with photon energies exceeding 4Ryd (54.4 eV) (Thuan 
|& Izotov|2005| ). Although a population of Wolf-Rayet stars has 


al.|2012t Sandstrom et al.|20r^ . Although the relatio n between 

COand H 2 has been the subject of much debate (see |Bolatto et 
|al.|2013| for a recent review), it is clear that at some low metal¬ 
licity, probably ^ 0.05 - 0.1 Z©, CO ceases to become a useful 
tracer of total H 2 mass. 

In this paper, we begin to explore the limits to which CO can 
trace H 2 in metal-poor star-forming ^laxies. We first describe 
the targets and ancillary data in Sect.J^ and the ^^CO(l-O) ob¬ 
servations and data analysis in Sect. ^ In Sect. we compare 
our results with previous work, and assess the dependence of CO 
luminosity on star-formation rate (SFR). Sect. [^describes con¬ 
straints on aco and Tdep, the molecular gas depletion time, i.e., 
the time required to consume the H 2 gas reservoir. We adopt a 
variation of the method proposed by [Saintonge et al.| ( |2011b| ), 
extending it to low metallicity in order to separate the effects of 
Tdep variations from variations of aco- In Sectl^ we analyze the 
H 2 content of our sample and compare it with H i for a compila¬ 
tion of galaxies taken from the literature. We discuss the impli¬ 
cations of our results in Sect. [7] and, in Sect. summarize our 
results and present our conclusions. 


2, The targets and the ancillary data 

The observing sample was defined mainly on the basis of de¬ 
tected H 2 emission, either in the ro-vibrational transitions at 
2/im, or the purely rotational transitions in the mid-infrared 
(e.g., Spitzer/IRS); however, we also targeted some galaxies 
that had already been observed in CO (see Sect. |3.4| ). Table 
lists the sample galaxies (ordered alphabetically by name), 
together with relevant observational parameters; Figs. and 
show thumbprint images. 


2.1. Individual source description 

The targets have several characteristics in common, in addition 
to their low metallicities. All of them, with the exception of 
NGC7077 and SBS 0335-052E, are sufficiently bright in the 
infrared (IR) to have been identified in the (faint) IRAS point- 
source catalogue. In fact, IIZw40 was identified as a 60-yum 


been detected in this galaxy ( [Zhang et al.|2007| ), these are insuffi¬ 
cient to power the high-excitation emission lines, thought instead 
to be due to fast radiative shocks ( [Thuan & Izotov|2005] ). |Hunt"^ 
ar] ( [2010[ ) found evidence for rotational H 2 emission, and despite 


the low metallicity of CGCG 007-025, also 7.7 pm and 11.3 pm 
emission from PAHs. Judgi ng from the relative ly high density of 
its ionized gas (-300cm"^, Izotov et al.|2007 ), CGCG007-025 


is an example of an “active” starburst, forming stars in compact 
dense regions ( [Hirashita & Hunt|2004[ ). 

★ IIZw40 (UGCA116) 

Although we did not observe this source in ^^CO(l-O), it is in¬ 
cluded in our CO survey of higher-/ lines (Hunt et al., in prep.); 
thus we consider it in our sample here. From its discovery, to- 
gether with I Zw 18, as an extragalactic Hit region ( [Sargent & 


Searle [T970j ), IIZw40 has been the subject of much scrutiny. 


Its two-tailed optical morphology and kinematics suggest that it 
is the result of a collision between tw o gas-rich dwarf galaxies 
which triggered the current starburst ( Baldwin et al.[[198^ [van 


[Zee et al7|[1998[ ). Near-infrared recombination lines show that 
more than half the 2 pm fiux comes from free-free emission, 
and warm dust contributes at least -25% to the 3-4 pm fluxes 
( [Joy & Lester 1988[ ). Emission from the ro-vibrational transi¬ 
tions of H 2 at 2pm was discovered by Vanzi et al. ( 1996[ ), and 
later found to be in a cloud distinct from the primary Hii region 
( Vanzi et al.|2008[ ). Visual extinction Ay is spatially v ariable in 
II Zw 40 and can be as high as ^10 mag in the nucleus ([Vanzi et 


[al.|2008[ ). This galaxy is undergoing a massive starburst in a re 
gion of -150 pc in diameter powered b y -5000-10000 OV stars 
( [Beck et ffi^[2002[ [Kepley et al.[|2014[ ). These massive stars are 
driving copious stellar winds which alter the observed kinemat- 
ics ([Bordalo et al.|2009[ ), and also produce high-ionization [Fev] 
( [Thuan & Izotov|2005[ ). Like CGCG 007-025, th e high density 


peaker source by [Vader et al.[ ( [1993] ). Except for N GC7077, all 
galaxies host conspicuous Wolf-Rayet populations ( [Schaerer ~et[ 
[al. 1 1 999}[Guseva et al.|2000[[Zhang et aT 2007[ ); they are all also 


derived from optical spectra (- 500 - 1400 cm~^, [Bordalo et al. 
2009[ ) and the small dimensions of the star-formation episode 
imply that it is an “active” starburst. 

★ Mrk996 

Mrk 996 is an unusual blue compact dwarf galaxy, with smooth 
early-type m orphology and a system of globular clusters ( Thuan 
1996[ ). Virtually all its star formation occurs in the cen¬ 
tral region of the galaxy (-315 pc in radius). Because of its 
high central density (- 10^cm ~^) and broad emission lines 


(-900 km s ^, Thuan et al. 1996), it was thought that there could 


be an accreting intermediate-mass black hole; however X-ray 
ob servati ons revealed no compact point source ( [Georgakakis et 


al 


2011). A high ionization [Oiv] line is detected in Mrk 996 
( [Thuan et al.[[2008[ ), that, like [Fev], needs -4Ryd to excite it. 


Thus, fast radiative shocks are also present in this galaxy, as its 


emission-line 
of ionization 


most metal-poor (12-rlog(0/H) = 7.7, 

-0.1 Z©) ^^CO(l-O) detection of our sample (see Sect. [^. 

Although it appears in an isolated galaxy catalogue ( [van Zee[ the sample of [ 
[2001[ ), it has a companion, a dwarf irregular galaxy at 8.3 kpc 
projected separation (84 km s“^ velocity separation) w hich may 

Zee[ 


properties are inconsistent with any other source 


Georgakakis et al. 2011). The lowest rotational 


transition of H 2 (-28 pm) is found in Mrk 996, the only one in 


iunt et al.] ( [2010[ ) in which it was detected. 


have triggered the current star-formation episode 


van 


200 1[). It is also found to be a member of a compact group by 


Focardi & KeE| ( [2002| ). CGCG 007-025 has th e bluest U -B 


color of any dwarf galaxy in the van Zee ( [2001 ) sample, with 


★ NGC 1140 (Mrk 1063) 

NGC 1140 is a barred dwarf irregular galaxy with peculiar mor¬ 
phology and a central starburst. The Ha luminosity of its cen¬ 
tral Hii complex is 40 times higher than that of the giant Hii 
region 30Doradus in the Large Magellanic Cloud. Within this 
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Fig. 1. Top row, from left to right in alphabetical order of name: SDSS composite (ugriz) images of CGCG 007-025, Haro 3 
(NGC3353), IIZw40, Mrk996. Bottom row, leftmost panel: NGC 1156 (image taken from DSS2, reprocessed); then from left 
to right SDSS images of NGC 7077, UM448, UM462. Images are oriented with North up. East left, and are 50" on a side, except 
for NGC 1156 which is ~10'. The color of IIZw40 is reddened by >2mag Ay from the Galaxy. NGC 1140 and SBS 0335-052E, 
both of which harbor SSCs and whose images are not from SDSS, are shown in Eig.|^ 


complex, roughly 50 0 pc in diameter, there are 6-7 blue super - 
star clusters (SSCs) ( [Hunter et al.||1994bt jde Grijs et al.||2004|), 
the most lum inous of these containing almost 6000 O stars ( jMolTj 
jet al.||2007| ). Although there is no galactic-scale outflow in this 
galaxy, powerful shocks from the massive SSCs are disrupting 
the IS M on scales up to 1-2 kpc radius ( jWestmoquette et ah] 
|2010| ). Kinematic and other evidence suggests that the formation 
of the SSCs and the larger-scale starburst have been triggered by 
an interaction or merger event with a low-luminosity late-type 
spiral about 1 Gyr ago ([Hunter et al.|1994a|bl [Westmoquette et[ 
[al.|2010[ ). 

★ NGC 1156 

A typical barred Magellanic irregular about 25% brighter in the 
B band than the Large Magellanic Cloud, NGC 1156 is the most 
quiescent galaxy in our sample. It has been described as the 
least disturbed galaxy in the Local Universe, with no signif¬ 
icant galaxy within 700 kpc ( [Karachentsev et slT. \1996\ . More 
recent H i observations revealed a single dwarf galaxy at a dis¬ 
tance of 80 kpc ( [Minchin et al.[MTU| ) which however is too 
faint to provide any measurable perturbation. The star forma¬ 
tion in NGC 1156 is spread out in Hii regions spanning the en¬ 
tire bar-rdisk, and the kinematics of the ionized (and atomic) gas 
and stars is consistent w ith circular rotation about a common 
axis ( [Hunter et al.|2002| . 

★ NGC 3353 (Haro 3, Mrk35) 

NGC 3353, the most metal-rich galaxy in our sample 
(12-rlog(0/H)~8.4, see Table [T]), is a dwarf irregular hosting a 
galaxy-wide starburst. However, high-resolution radio images 


show that most of the starburst activity is concentrate d within 
a single knot to the northwest of the galaxy nucleus (Johnson 
[et al. [2004] ). The star clusters in this knot are quite young (1- 
4Myr) and massive (~ 10^ M©) and reside in a region with 
Ay~2-8 mag, implying that much of the current star formation 
is currently enshrouded in dust ( [Johnson et al.[[20()4| ). In fact, 
optical measures of the extinction are smaller than those mea¬ 
sured with the radio and n ear-infrared observations ( [Steel et 
1996 ; Cairos et al. 2007[ ). Detailed photoionization models of 


the infrared fine-structure lines suggest that there are at least two 
emission zones in this galaxy: one with low extinction, visible 
in the optical, and another, optically invisible with much higher 
extinction ( [Hunt et ar|[2006[ ). Rotational H 2 has been detected 
in NGC 3353, and the detection of the high-excitation [Oiv] line 
(^4 Ryd) implies that radiative shocks are also present ( [Hunt et 
al.|2010 ). Like other galaxies in this sample, kinematic evidence 


suggests that a merger is responsible for triggering the starburst 
in NGC 3353 ( [Cairos et al.|2'007l ). 

★ NGC 7077 (MrkOOO) 

NGC 7077 is classified as an early-type galaxy, E/SO, but with 
peculiar morphology, and relatively blue colors ( [Gil de Paz et 
[ar][2003[ Micheva et al.[[2013] ). It is an isolated galaxy, part of 
the [van Zee[ ( [2000[ ) sample, with red colors in its outer regions 
(see also Lig.[lJ). These red colors are in contrast with the strong 
(blue) nuclear s tarburst which has an intensity comparable to that 
of 30Doradus ( van Zee|2000[ ). Perhaps the most salient feature 
of NGC 7077 is its extended ultraviolet disk, suggesting that the 
disk of this galaxy is in an active growth phase ([Moffett et al. 

[MT^ . 
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Table 1. Observing targets 


Name 

Other 

Redshift^ 

Systemic 

Distance's 

Distance's 

12+ 

Dimensions^s 

SFRs^ 

Logs^ 

Logs^ 

Log(Total)® 


name 


velocity‘s 

(Mpc) 

method 

logiO/Hf 

(arcsec^) 

(Moyr-i) 

Mstar 

Mm 

Mhi 




(kms“^) 






(Mo) 

(Mo) 

(Mo) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

CGCG 007-025 


0.0048 

1449 

24.5 

CMB 

l.lAr 

27.0 X 16.2 

0.23 

8.36 

8.90 

8.90 

IIZw40 

UGCA116 

0.0026 

789 

11.7 

CMB 

8.12 

33.6 X 13.2 

1.75 

8.25 

7.54 

8.77 

Mrk996 


0.0054 

1622 

18.1 

CMB 

8.00 

36x30 

0.16 

9.18 

8.00 

8.00 

NGC 1140 

Mrk 1063 

0.0050 

1501 

19.7 

TF 

8.18 

102 X 54 

0.82 

9.76 

8.56 

9.51 

NGC 1156 


0.0013 

375 

8.1 

Stars 

8.23 

198 X 150 

0.23 

9.55 

8.66 

9.04 

NGC 3353 

Haro 3, 

Mrk 35 

0.0031 

944 

18.1 

TF 

8.37 

72.0 X 49.8 

1.28 

9.54 

8.87 

8.87 

NGC 7077 

Mrk 900 

0.0038 

1152 

17.2 

TF 

8.03 

48 x42 

0.12 

8.98 

8.15 

8.15 

SBS 0335-052E 


0.0135 

4043 

53.6 

CMB 

7.25 

13.8 X 12.0 

1.30 

8.51 

7.81 

8.61 

UM448 

Arp 161, 
UGC 6665 

0.0186 

5564 

81.2 

CMB 

8.04 

24x24 

11.11 

10.91 

9.84 

9.84 

UM462 

Mrk 1307 
UGC 6850 

0.0035 

1057 

19.5 

CMB 

7.97 

36x30 

0.29 

8.56 

8.30 

8.72 


NASA/IPAC Extragalactic Database (NED, http://nedwww.ipac.caltech.edu/): CMB corresponds to Cosmic Microwave Background and TE to 

Tully-Eisher; velocities are heliocentric; dimensions are optical major and minor diameters. __ 

Direct-temperature metallicities for CGCG 007-025, NGC 3353, UM 448, UM 462, SBS 0335-052E from|lzotov et al.|(|2007|); II Zw 40 (|Thuan 
& Izotov|2005} ; Mrk996 jThuan et al.|l996) ; NGC 1140 ( jlzotov & Thuan|2004t ; NGC 1156 jVigroux et al.|1987t ; NGC 7077 j van Zee & 


Haynes |2006| ). 


These quantities are newly derived, as described in Sect. |2.3| 

When these values differ from those in Col. (12), they are calculated from the mean Hi surface density integrated over the optical size of the 

galaxy. _ __ 

The Hi measurements are taken from |van Zee|{|2001| CGCG007-025, NGC7077), |van Zee et al.|(|1998| IIZw40, UM462),|Thuan et al.| 
(fT999l Mrk996),|Hunter et al.|(|1994a| NGC 1140), Swaters et al.|([2^ NGC 1156), |Hunter et al.|(|1982| NGC 3353), |Ekta et al.|(|2009| 
SBS 0335-052E), lPaturel et al.lp003r UM448). 


★ SBS0335-052E 

Like, IIZw40, we did not observe this source in ^^CO(l-O), 
but it is included in our CO survey of higher-/ lines (Hunt et 
al., in prep.); thus we consider it here. This is the most metal- 
poor galaxy of our sample (12-rlog(Q/H )~7.2), with a neigh - 
bor at ~22kpc distance (SBS 0335-052W, |Pustilnik et al.|2001| ). 


This companion, embedded in a common H i cloud, is the most 
metal-poor star-forming galaxy known in the Local Universe 


(12-rlog(0/H)~7.1, Lipovetsky et al.[[1999l 

Izotov et al.[[2005 ). 

High-resolution optical and near-infrared i 
the star formation activity in SBS 0335-052 
massive SSCs, with the two youngest ones ( 
prising -5000 07.5 V stars (Reines et al.|2C 
tremely low abundance, SBS0335-052E is 
tremelv high extinction (Av^ 15 mag, [Plant 
Houck et al.||2004), and ^ 10"^ M© of dust si 
two youngest SSCs ([Hunt et al.|2001 2014J. 
present in this galaxy, with strong emission 
transitions at 2yum ([Vanzi et al.|2000 [2011). 

maging shows that 
H is confined to six 
:^3 Myr) each com- 
)08|). Despite its ex- 
known to have ex- 
e & Sauvage[[2002l 

jrrounding only the 
Warm H 2 is clearly 
from ro-vibrational 
, The cm-regime ra¬ 
ng ionized-gas den- 

dio spectrum is highly self-absorbed, imr 

)lyi 

sities of ^3000cm ^ (Hunt et al. 

2004 

Johnson et al. 2009). 


SBS 0335-052E is the proto-typical “active” starburst with its 
star-formation activity occ urring in compact (unres olved at HST 
resolution), dense regions ( [Hirashita & Hunt|2004| ). 

★ UM 448 (Arp 161, UGC 6665) 

UM448 is a merger remnant with residual tidal tails, and so IR 
luminous (Ljr = 6.3 x 10^^ L©) that it would almost be clas¬ 


sified as a Luminous IR Galaxy (LIRG, [Dopita et al. 2002). 
Despite its low metal abundance (12-rlog(Q/H)~ 8.0, 0.2 Z©, see 
Table [^, this galaxy shows clear PAH emission ( Galliano et al. 


2008| ). This is perhaps because of a metallicity gradient which 
reaches a peak value of 12-i-log(0/H)~8.4 toward the northwest, 
a region of lower temperature ( [James et~ar]|2013| ). There is ev¬ 
idence for age gradients in the stellar populations and complex 
ionized emission-line kinematics, both of which probably reflect 
the merger nature of UM 448 ( [James et al.|2013| ). 

★ UM462 (Mrk1307, UGC 6850) 

Like many of the galaxies in our sample, UM462 is possi¬ 
bly tidally interacting with a relatively close neighbor, UM461 
( [Taylor et al.[ 1995[ ), although [van Zee et al.[ ( [199^ concluded 
that the Hi kinematics are inconsistent with this. Most of the 
star formation in UM 462 occurs in two main knots of star for¬ 
mation, although high-resolution observations res olve at least 
four starbursting regions with disrupted kinematics ( Vanzi et al. 
2002[[James et al.|201Q[). The Wolf - Rayet content of this galaxy 


has been disputed by [James et ah] ( [2010[ ) who do not fi nd evi 
dence to su pport a significant Wolf-Rayet classification. Vanzi 
et al. ( 20Q2[ ) found ro-vibrational H 2 emission in this galaxy, but 


no evidence for [Feii], implying that the star-formation episodes 
must be younger than ~10Myr. The knots of star formation are 
so compact and luminous as to be classified as SSCs ([Vanzi 
[ 200 ^ . 

2.2. An additional source, (Zw 18 

Although IZw 18 was not included in our sample, it was ob- 
served, although n ot detected, in the ^^CO(l-O) transition by 
Leroy et ^ ( |2007| ). Because of its extreme low metallicity, sim¬ 


ilar to SBS 0335 -052E, we have inco rporated this galaxy in our 
analysis (see also[Schruba et al.[2012[). We have adopted the SER 
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from |Hunt et al.| ( |2005] ), and calculated the stellar mass, Mstar, as 
described below (see also |Fumagalli et al.|2010| ). 

Like IIZw40, IZwl 8 was discovered as an extragalac- 
tic Hii region by |Sargent & Searle| (|1970|). Together with 
the SBS0335-052E and SBS0035-052W pair, it is the most 
metal-poor star-forming galaxy in the Local Universe with 
12-rlog(0/H) = 7.17-7.22 ( Thuan & Izotov|[2()05| ). Most, if not 
all, of the current star formation occurs in two main clusters, 
although the “C component” (or Zwicky’s flare) is also form¬ 
ing stars ( [Izotov et al.|2001| ). However, the two clusters are not 
as extreme as those in SBS0335-052E (and other galaxies in 
our sample), being somew hat less luminous than 30Doradus 
( [Hunter & Thronson||1995| ). Like SBS 0335-052E, around the 
main star clusters there is a complex structure of filaments 
and arcs from ionized gas, and the complex Hi morphology 
and kinematics suggest that the current star-formation ep isode 
was triggered by an interaction with the C component ( [Lelli[ 
[et al.[[201^ . Unlike SBS 0335-052E, however, there is no ev¬ 
idence for a dense gas component in IZw 18; the radio contin¬ 
uum spectrum follows the trend expected for optically thin free- 
free emission ( [Hunt et al.|2005[ ) with an ionized-gas density of 
:^ 100 cm“^, similar to that found from optical spectra (Izotov et 


al.|2001[ ). Despite their similarly low metal abundance, the dust 


masses, dust-to-gas, and dust-to-stars mass ratios in IZw 18 are 
more than 100 times lower than in SBS 0335-052E ( [Hunt et ak] 
[2014[ ), implying that the density of the cool gas in the ISM may 
play an important role in dust formation ( [Schneider et al.|20T5 ). 



Fig. 2. NGC1140 and SBS 0335-052E taken from composite 
HST images; both galaxies h ost super-star clusters ( [Hunter et 


al.|1994bt[Reines et al. 12008] ). The image of NGC 1140 is ~90 
on a side, and SBS 0335-052E ~5", and orientation is North up. 
East left. 


2.3. SFRs and stellar masses 


We ha v e calcu lated the SERs for our targets according to [Calzettr 
et al.[ ( [2010[ ) using Ha and 24yum luminosities (their Eqn. 


17). This approach obviates the need for uncertain extinction 
corrections, and has the added advantage of accounting for 
the general behavior of dust opacity in galaxies. Many of the 
Ha total fluxes come from [Gil de Paz et al.[ ( [2003 [ Haro 3, 
IIZw40, Mrk996, NGC 7077), but others come from [van Zee[ 
( [20001 CGC G 007-025), [Kennicutt et al.[ ( [200^ NGC 1140, 


NGC 1 156), [bopita et al.[ ( [2002| UM44'8^), and 


James et al. 


(2010 UM 462). The 24 yum fluxes are taken from Engelbracht et 


[ar]( [2008l ), [Tbmi et al.[ ( [200^ NGC 7077), or Hunt et al. (2015, in 
prep.). Eor SBS 0335-052E and IZw 18, we take the SER mea- 
sured by high-frequency radio free-free emission ([Johnson et al. 
[200^[Huntetal.|2005[). _ 

Eollowing [Grossi efST ( [2015 ), stella r masses, Mstar, w ere es¬ 
timated according to the recipe given by Wen et al.] ( [2013[ ) based 
on WISE W1 (3.4yum) luminosities. This approach exploits the 
relatively constant mass-to-light ratios of stellar populations, in¬ 
dependently of metallicity and age, at near-infrared wavelengths 
( [Norris et al.[[2014t [McGaugh & Schombert[[2014| ). However, 
when W1 photometry was not available, we used IRAC 3.6yum 
photometry instead. In fact, the two bands are very similar; us¬ 
ing data from [Brown et al.[ ( [2014[ ), [Grossi et~ar] ( [2015[ ) And for 
spirals a mean flux ratio F 3 . 4 /F 3.6 = 1.02 + 0.035. Incl uding 
also the data for dwarf irregulars from [Brown et al. ( [2014[ ) gives 
a mean flux ratio F 3 . 4 /F 3.6 = 0.98 + 0.061. Thus, we conclude 
that the ratio of the W1 and IRAC 3.6yum bands is unity, with 
5-6% scatter for galaxies like our targets. 

We used the Hii-galaxy formulation by [Wen et aL| ( [2013[ ), 
rather than what they found for their full sample; the Hii galaxies 
have the lowest mass-to-light ratios in their compilation, corre- 
sponding roughly to the bluest regions of the galaxies studied by 
[Zibetti et"ar] ( [2009[ ). To better take into account the weak trends 
with abundance found by [Wen et al.[ ( [2013] ), we also applied an 
approximate correction f or low metallicity (by multiplying the 
mass-to-light ratio by 0 . 8 : f^n et al.|20T3 see their Eig. 17). 


However , before applying the Hii-galaxy relation by [Wen 
[et al.[ ( [2Q13[ ), we first subtracted nebular emission and emis¬ 
sion from hot dust. Such contamination can be very important 
in the near-infrared and can cont ribute 50% or more to the ob- 
served flux at these wavelength s ([Hunt et al.|2Q(3T 20Q2[ Smith 


& Hancock[[2009t [Hunt et ^[26l2[ ). The ionized gas contin¬ 


uum contribution to the 3.4-3.6yum flux was estimated from the 


SER using the emission coefficients from Osterbrock & Eerland 
( [2006[ ). The subtraction of hot dust is based on the assump¬ 
tion that i7-band emission is entirely stellar. Because i7-band 
photometry is available for most of our sample, we used the 
data from Brown et al.[ ( [2014[ ) to estimate the maximum pos¬ 
sible IRAC 3.6yum/i7-band ratio in galaxies similar to our tar¬ 
gets; 95% of the spiral/dwarf irregular galaxies have a flux ratio 
<2.4. This corresponds to a (Vega-based) [H-3.6] color of ~ 0.8, 
consistent with what is found for the pure stellar component in 
star-forming galaxies ( [Hunt et al.[2002| ). After subtraction of the 
nebular component, any excess over this ratio was attributed to 
hot dust and subtracted. We found that hot dust can contribute 
between 10% and 50% in the 3.4-3. 6 yum bands (most signifi¬ 
cantly for IIZw40, SBS 0335-052E). 

We compared the stellar masses obtained with the formu- 


lation of 

Wen et al. (2013[) to those calculated according to 

Lee et al. 

d2006 

) based on IRAC 4.5 jum luminosities (used by 
), and And they are consistent to within 0.2 dex. 

Hunt et al.|2012 


Moreover, for the galaxies having both W1 and IRAC data, the 
stellar masses obtained from 3.6yum luminosities are within ~5% 
of those from W1 as expected. The Kroupa-like broken power- 
law initial mass function (IME) used by [Calzetti et al. ( 2010| ) and 
consequently for our calculation of SERs is consistent with the 
IME for Mstar by [Wen etal.| ( |MT^ . 


2.4. HI measurements 

All the sample galaxies have been observed in H i, half of them 
with interferometric maps. Eor the galaxies with resolved obser- 
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vations, for comparison with CO, we have attempted to compute 
the Hi mass, Mhi, by integrating the mean Hi surface density 
within only the optical radius. Because virtually all metal-poor 
dwarf galaxies show H i emission well outside their optical ra¬ 
dius (e^g^jHunteFeF^D[T994ayya^^Zee[£t^D[^^ |van Zee 


|2000t |Ekta et al.||2009||Lelli et al.||2012||2014| r these values of 
Mhi are typically a factor of 2-10 lower than the total Hi mass, 
integrated over the entire Hi extent. Tablereports both values 
of Mhi when resolved H i maps are available. 


2.5. Comparison samples 

To expand the parameter space for our analysis, w e also con- 
sidered another dwarf sample with CO detections ( [Schruba 
aL]|2012| ), as well as no rmal star-forming galaxies, LIRGs, and 


Ultra-LIRG s (ULIRGs, Gao & Solomon| 2004 |Gracia-Carpio| 


|et al.|[200^ |Garcia-Burillo et al.||2012 ). For the latter galax¬ 
ies, we converted their published IR luminosity to SFR ac¬ 


cording to Kennicutt] ( 1998 ), after applying the correction to 


Kroupa (|26oi 


Schruba et al. 


IMF for consistency with our SFRs and those of 
( [2012| ). For the Schruba et al. CO(2-l) observa¬ 
tions, we applied the factor they recommend to convert ^^CO(2- 
1) to ^^CQ(l-O): Ico 2 -illcoi-o = 0.1. Stellar masses for the 


Gao & Solomon|(|2004[) and|Garcia-Burillo et al.| ( |2012[ ) samples 
were taken fro m|Skibba et al.| ( |2011| ), and for the galaxies from 
Schruba et al.|(|2Q12|), masse s were adopted either from [Hunt et| 
al.| ( |2012| ) or |Lee et ST| ( |2006| ). These latter estimates are based on 
IRAC photometry similar to the method used here, as described 
in Sect. 12.31 

We incorporated also the galaxies from the CO(2-l) APFX 
Low-redshift Legacy Survey of MOlecular Gas (ALLSMOG, 
Both well et al.| 2014| ). This sample is essentially an exten¬ 
sion to lower stellar masses (8. 5 < log(Mstar/M0) < 10) o f the 
COLDGASS gal axies studied byjSaintonge et al.| ( |2011a|b| ). For 
consistency with Schruba et al. ( 2Q12| , to convert CO(2-1) 
to ^^CO(l-O), we again u se /co2-i/-^cqi-o = 0-7 , rather than 
^co2-i/^coi-o = 1 as in [Bothwell et al.| ( |2014| ). Because of 
the -0.4 dex offset between the ALLSMOG metallicity cali¬ 
bration an d the mainly direct-temperaturj^ determinat ions used 
here (e.g.. Hunt et al.|2Ql 2} [Andrews & Martini|2Q13| ), the oxy¬ 
gen metallicities given by Bothwell et al. | ( |20 1 4| ) wer e converted 
to the[Pettini (^Pag^ ( |2004| ) calibration according to Kewley ^ 


Ellison 


2008| |^ Finally, to extend the metallicity range to ex 


trem ely metal-poor gala xies (XMPs), we incorporate the sample 
from |Filho et ST] ( [2013| ) which, however, have no CO measure¬ 
ments; we include them to better investigate gas content at low 

metallicities. _ 

Both the ALLSMOG and the |Filho et 3T] ( |2013| ) galaxies 
have stellar masses and SFRs taken from the MPA-JHU cat¬ 
alogue ( [Kauffmann et al.||2Q03t [Brinchmann et al.j|2Q04| ), the 
same basis for calibration used by | Wen et al. ( |2013| . However, 
the automatic algorithm used by the SDSS does not always pro¬ 
vide the correct global photometry for galaxies whose surface 
brightness is not uniformly distributed (e.g., in clumps or knots). 
This should not be a problem for the ALLSMOG galaxies, which 
tend to be extended spirals, but cou ld be a problem for the com¬ 
parison with the XMP sample by |Filho et ST] ( [2013| ). In fact, 

^ The electron temperature derived from the ratio of the [Oiii] 24363 
auroral line to lower-excitation lines ([Oiii] 24959,5007) is converted 
to a metallicity after correc ting for unseen phases of ionization (e.g., 

Osterbrock & Ferland|2006|). _ _ 

^ These metallicites are shown in Figs.[5|[^|8||^[ll|[T^ and 14 


Filho et ^ ( 2013| ) found that comparing mass-to-light ratios 


based on optical colors gave as much as an order of magnitude 
Mstar excess over the SDSS estimates; the discrepancies tend to 
be largest for galaxies with either cometary or two-/multiple- 
knot morphology. To check for systematic differences between 
the near-infrared method used here (IRAC, WISE) and in |Hunt| 
et al.| ( |2012| ) (see also [Madden et al.||2Q13] ) and the optically- 


determined (MPA-JHU) SDSS values, we searched the MPA- 
JHU database and found 17 low-metallicity galaxies for which 
we had calculated stellar masses in [Hunt et ^ ( [2012 ) using 
4.5 pm luminosities. Although the contributions of the nebular 
continuum and hot dust emission have been subtracted, the val¬ 
ues of Mstar using IRAC/WISE photometry are, on average, 17 
times larger than those with SDSS with a scatter of -0.4 dex. We 
investigated the morphology of the galaxies in this comparison 
and find that they tend to be very clumpy, with either multiple 
knots or cometary features. Comparing the SDSS photometry 
with the broad-band JHK photometry for these galaxies gives 
NIR-optical colors that are extremely red, unrealistically so for 
these blue galaxies; the implication is that the optical fluxes are 
being underestimated. We thus ascribe this difference to the diffi¬ 
culty of the SDSS automated photometric algorithms when con¬ 
fronted with multiple clumps. As a consequence, in our compar¬ 
isons with the XMP sample, these uncertainties should be kept 
in mind. 


3. The CO observations 

We have observed ^^CO(l-O) in our sample galaxies over a 
three-year period from 2008 to 2010 with the IRAM 30-m tele¬ 
scope (Pico Veleta, Spain). 


3.1. IRAM observations 

In September, 2009, and September, 2010, we observed 
six galaxie^ in ^^CO(l-O) with the Eight Mixer Receiver 
(EMIR) using the Wideband Line Multiple Autocorrelator 
(WILMA) backend. Observations of 1^00(1-0) for NGC 1140 
and NGC 1156 were performed with the AB receivers in 
November, 2008, using the 1 and 4 MHz backends at 3 and 
1.3 mm (88-115 GHz, 219-229 GHz, respectively). The precip- 
itable water vapor was ^4 mm during the 2008 run, giving sys¬ 
tem temperatures of Tsys of :^370K at 3 mm. In 2009, the con¬ 
ditions were worse (^11 mm), but the EMIR receivers had com¬ 
parable Tsys, -335 K at 3 mm; in 2010, our most extensive run, 
Tsys varied from -230 K to ^400 K, with precipitable water va¬ 
por ranging from - 4 mm to 10 mm. In all runs, we used wob¬ 
bler switching with a throw of 90". Intensity calibration was per¬ 
formed with two absorbers at different temperatures. 

Pointing was checked at regular intervals (every 1.5 hr) on 
nearby planets or bright quasars (e.g., 0234-r285), and focus 
tests were carried out every 4 hrs during the night and every 3 hrs 
during the day. Pointing errors never exceeded -3", and were 
usually ^2". Given the beam size (half-power beam width) of 
-22 ", this was more than sufficient. 


IIZw40 and SBS 0335-052E were observed only in the higher-J 
lines with APEX. 
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Fig. 3. ^^CO(l-O) average spectra for the observed galaxies. Top row, 
from left to right in alphabetical order of name: CGCG 007-025, Mrk996, 
and NGC 1140. Second row: NGC 1156, NGC3353 (Haro 3) two pointings 
([0,0], [0,+10"]. Third row: NGC 3353 third pointing [-6", -9"]), NGC 7077 
(Mrk 900), UM 448. Bottom row: UM 462. The baselines are shown as a hor¬ 
izontal solid lin e, to gether with the multiple-component Gaussian fits as de¬ 
scribed in Sect. |3.2[ The vertical axes are in T^ units, while in Table the 
units are Tnib, converted from T^ as described in the text. Velocity channels 
are as reported in Table 


3.2. Data reduction 

We adopted the GILDAS/CLASS data reduction package to ob¬ 
tain mean spectra (http://www.iram.fr/IRAMFR/GILDAS). To 
remove the continuum from each spectrum, a polynomial base¬ 
line was fitted to the line-free regions of each scan (defined 
a priori) and subtracted; the scans were thereafter smoothed 
and averaged, and a constant baseline subtracted. We measured 
the peak intensities, central velocities, full width half-maximum 
(FWHM) and velocity-integrated fluxes of the detected lines by 
fitting Gaussian profiles to the data. In most cases, there is sig¬ 
nificant velocity structure in the emission, so we fit the line pro¬ 
files to multiple Gaussians. In these cases we use the sum of the 
integrated line intensities from the multiple Gaussian fits in the 
analysis (see Table |^. 

Antenna temperatures (T^ have been converted to main- 
beam temperatures (Lmb) by dividing the antenna temperatures 
hy T] = ^eff/Feff, where ^eff (=0.78) and Fgff (=0.94) are 
the beam and forward hemisphere efficiencies, respectively (i.e., 
Tmb = Ty77). To convert the measured ^^CO(l-O) brightness 
temperatures, [Lmb (K)] to fluxes [S (Jy)], we used the standard 
IRAM 30-m conversion factor of 3.906 JyK“^ The ^^CO(l-O) 
line profiles are shown in Fig. together with the Gaussian fits 
obtained using the GILDAS/CLASS package as described above 
and reported in Table Although a few of the individual veloc¬ 
ity components are too narrow to be significant (e.g., NGC 1140, 


NGC 1156), we have checked that the total intensities are robust 
to the details of the fitting. 


3.3. Aperture corrections 


The IRAM 30 m has a beam with an FWHM of ~22 " at 3 mm. 
Consequently, the largest galaxies in our sample (NGC 1140, 
NGC 1156, NGC 3353, see Table may be missing flux. We 
have attempted to corre ct for this by adopting the approach of 
|Saintonge et al.| ( |2011a| ) who calibrated aperture corrections us¬ 
ing the set of nearby g alaxies observed by|Kuno et al.| ( |20()7| ). 
We applied Eqn. (2) by [Saintonge et al. ( |2()lla ) based on opti¬ 
cal size, but used the geometric means of the major and minor 
axes given in Table These corrections amount to an increase 
in line intensities by a factor of 1.64 for NGC 1140 and 3.23 for 
NGC 1156. 

For NGC 3353 (Haro 3), we have additional pointings spaced 
by roughly 0.5 beams, hence not independent. Where necessary, 
Saintonge et al.| ( [2011a| ) used spacings of 0.75 beams along the 


major axis, and found a mean ratio between the offset flux and 
the central flux of 0.33. The corresponding ratio for the [-6", 
-9"] pointing of NGC 3353 (along the southern major axis) is 
0.53, larger than the offs et corresponding to the wider spacing by 
[Saintonge et al.|p011a| ), but consistent with what could be ex¬ 
pected for the closer sp acing of our offset. Hence , for NGC 3353, 
we adopted Eqn. (3) of [Saintonge et al. ( |2011a[ ) with their mean 
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Table 2. ^^CO(l-O) emission line parameters 


Name 

Telescope time 
(min) 

Component 

Channel 

(kms"^) 

Tr„s(mb)‘‘ 

(K) 

Offset 

(kms"^) 

FWHM 

(kms"^) 

Integrated intensity"^ 
(Kkms"^) 

CGCG007-025 

524 

Total 

10.5 

0.0018 

-39 

115 

0.286 (0.095) 



Total 

15.7 

0.0016 

-39 

105 

0.285 (0.102) 

IIZw40" 

- 







Mrk996 

573 

1 

5.2 

0.0019 

-1 

28 

0.18(0.03) 



2 

5.2 


39 

10 

0.05 (0.02) 



Total 

5.2 


0 

- 

0.23 (0.04) 

NGC1140 

831 

1 

5.2 

0.0016 

-17 

26 

0.32 (0.05) 



2 

5.2 


3 

8 

0.07 (0.04) 



3 

5.2 


30 

38 

0.18 (0.05) 



4 

5.2 


64 

7 

0.04 (0.02) 



Total 

5.2 



- 

0.61 (0.08) 



Total^ 

5.2 


0 

- 

1.00 (0.13) 

NGC1156 

47 

1 

10.4 

0.0076 

-84 

13 

0.36 (0.09) 



2 

10.4 


-35 

100 

0.89 (0.09) 



3 

10.4 


44 

42 

0.48 (0.09) 



Total 

10.4 


0 

- 

1.25 (0.16) 



Total^ 

10.4 

0.0076 

0 

- 

4.04 (0.52) 

NGC3353 [0",0"] 

79 

1 

5.2 

0.0056 

-31 

60 

1.83 (0.17) 



2 

5.2 


21 

16 

0.34 (0.11) 



Total 

5.2 


0 

- 

2.17(0.20) 



Total^ 

5.2 


0 

- 

3.07 (0.29) 

NGC3353 [0",-hl0"] 

77 

1 

5.2 

0.0045 

-35 

42 

1.38 (0.12) 



2 

5.2 


19 

22 

0.43 (0.12) 



Total 

5.2 


0 

- 

1.81 (0.17) 

NGC3353 [-6",-9"] 

80 

1 

5.2 

0.0046 

-17 

71 

0.84 (0.20) 



2 

5.2 


25 

36 

0.62 (0.18) 



3 

5.2 


98 

31 

0.30 (0.08) 



Total 

5.2 


0 

- 

1.46 (0.27) 

NGC7077 

585 

1 

5.2 

0.0017 

-5 

20 

0.25 (0.03) 



2 

5.2 


35 

30 

0.20 (0.04) 



Total 

5.2 


0 

- 

0.45 (0.05) 

SBS 0335-052E" 

- 







UM448 

454 

1 

5.3 

0.0017 

-115 

20 

0.11 (0.04) 



2 

5.3 


-58 

63 

0.42 (0.16) 



3 

5.3 


1 

44 

0.52 (0.14) 



Total 

5.3 


0 

- 

1.04 (0.22) 

UM462 

119 

Total 

10.4 

0.0034 

-9 

121 

0.824 (0.173) 



Total 

15.7 

0.0032 

-8 

116 

0.749 (0.207) 


This and all temperature units in this table are main-beam, Tmb- 

Values in parentheses are the uncertainties. The total line flux is the sum of the individual components, and the total uncertainty is calculated by 

adding in quadrature the individual uncertainties. _ _ 

As discussed in the text, we did not observe this galaxy in ^^CO(l-O), but rather relied on the data reported by Taylor et al. ( |1998| ). 

With aperture corrections as described in the text. 

We did not observe this galaxy in ^^CO(l-O), but rather adopted the ^^CO(3-2) data reported by Hunt et al. ( 2014 ), converted to CO(l-O) 
assuming a CO(3-2)/CO(l-0) flux ratio of 0.6. 


ratio between offset and central intensities (0.33); this correc¬ 
tion gives a multiplicative factor of 1.41. Tablealso reports the 


^^CO(l-O) intensities corresponding to the aperture corrections 
for these galaxies which we will adopt throughout the paper. 
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3.4. Comparison with previous CO observations 

Because some of the galaxies in our sample have pr evious CO 
observations, here we compare our results with those. [Sage et al 


( [T99^ observed NGC3353, NGC7077, UM448, and UM462, 
and securely detected only the former two. The results for 
NGC3353, obtained with the IRAM 30-m, are in good agree¬ 
ment with our central pointing (they find an integrated CO in¬ 
tensity of 1.98+0.28 Kkms“^), although we find ~70% of their 
CO emission for NGC7077 (they find 0.68+0.1 Kkms“^). Our 
detections of UM448 and UM462 are consistent with their 
marginal detections obtained with the NRAO 12-m. 

UM462 was also observed by [Brinks & Taylorj \\99l\ 
and [Gondhalekar et al.] ( [1998] ), neither of whom detected the 
galaxy; their upper limits (<0.28Kkms“^ and <0.5Kkms“\ 
respectively) are somewhat discrepant with our detecti on (0.7 
0.8 Kkms“^), as well as with the marginal detection of Sage et 
Z1([1992| ~0.7Kkms-i) 



z 1 1 

1 iiiii[ — r 

1 1 iini| — 1 

1 1 1 llll| - 1 1 1 IIHIj— 

1 1 1 lllll| -1 1 1 lllll| 

□ r 

L 10* 
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. •"vHwd"" 


iiiii| 


([1992| ~0.7Kkms- __ 

NGC 1140 was observed previously by [Albrecht et al.[ ( f2004| ) 
and NGC 1156 by [Leroy et al.[ ( [2005] ). The two previous obser¬ 
vations are reasonably consistent with our results. 


4. CO luminosity and SFR 
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Our new IRAM observations have successfully detected 
^^CO(l-O) down to a metallicity of 12+log(0/H) = 7.74 
(CGCG 007-025, signal-to-noise ratio 3.0, see Table [^, which 
is the lowest metallicity at which CO has been (tentatively) de- 


tected (see WLM: 12+log(0/H) = 7.8, [Elmegreen et al.j[2()13| ). 
To our knowledge. Table [^ gives the largest number of ^"^CO(l- 
0) detections ever reported for 12+log(0/H)^8 (0.27©) outside 
the Local Group. 

In the following, we compare the integ rated CO luminosity, 
L'with SFR as calculated in Sect. 


2.3 


L^ was derived ac¬ 


cording to the formulation of [Solomon & Vanden Bout ( [2005| ), 
using the appropriate Jy-to-K conversion (see Sect. [3 for the 
IRAM value). The result is shown in Fig.[^ 

Not surprisingly, most of the metal-rich star-forming galax¬ 
ies in the comparison samples follow the linear correlation be 



tion but offset to lower L^q by a factor of ~30: 

log(SFR) = log(L^o) - (6.92 + 0.07). 

The best-fitting slop e obtained with a robust fitting technique 
( Li|1985[[Fox|1997[ ) i s unity within the errors ( 1.05+0.15); thus, 
for consistency with [Gao & Solomon[ ( 200^ , we fixed a unit 
slope and derived the corresponding intercept as give n above and 
sho wn in Fig.[^ Although we did not use the data of [Schruba et[ 
[ak ( [2012| ) to establish the best-fit regression, most of those galax¬ 
ies also follow the trend for our sample (see Fig.[^. The excep¬ 
tions are the two most metal-poor galaxies, SBS 0335-052E and 
IZwl8 . 

Judging from our sample alone, the offset between SFRs 
and high L^q and low L^q is apparently independent of metal¬ 
licity, at least to the lowest oxygen abundance of our sample, 
12+log(0/H)~7.7. However, CGCG 007-025 at this metallicity 
is only tentatively detected, and our other detected galaxies have 


Fig. 4. SFR vs. L^q for our sample a nd a dditional galaxies from 
the literature. As discussed in Sect. |2.2| we include IZw 18 as 
one of our sampl e galaxies (see also Schruba et al.[[2012| ). The 
open triangles for Schruba et al. ( [2012 ) correspond to their com¬ 
pilation of literature data and the filled ones to their metal-poor 
dwarf observations. CO non-detections are shown by leftward 
arrows. The rightmost regression line corresponds to the linear 
correlation found by [Gao & Solomon] ( 2004| ), and the leftmost 
one to our new determination, also compatible with a slope of 
unity (see text). The horizontal dotted line gives the IR l uminos- 
ity limit above which the linear LIR-L^q correlation of Gao & 
Solomon ( 2004| ) no longer holds. 


a relatively narrow spread in metallicities, 12+log(0/H)~8. If we 
include in the regression the [Schruba et al.[ ( [2012b compilation 
and their five ^^CO(2-l) detections (converted to '^CO(l-O) as 
described above), we find a significant correlation between the 
SFR/L^q ratio and 12+log(0/H), as shown in Fig. [^ (Pearson 
parametric correlation coefficient, -0.88, >99.99TO signifi¬ 
cance level). Two galaxies are discrepant with this metallicity 
trend, CGCG 007-025 and SBS 0335-052E. SBS0335-052E 
was not directly observed in ^^CO(l-O), so there may be some 
problems with the multiplicative factor relating CO(3-2) to 
CO(l-O) (see [Hunt et aklpT^ . CGCG007-025, on the other 
hand, is directly observed in ^^CO(l-O) but the detection is 
only tentative. More CO(l-O) observations at metallicities be¬ 
low 12+log(0/H) = 8 are needed to verify this inconsistency, es¬ 
pecially because of the expectat ion that for Z/Zq^0 .1, CO may 
cease to be a viable tracer of H2 ( [Bolatto et al.|2013| ). We discuss 
the implications of the relation between SFR/L^q and abundance 
in the following section. 
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18+log(0/H) 


Fig. 5. SFR/L' vs. 12+log(0/H) for our sample and addi- 


tional galaxies from the literature. As discussed i n Sect. |2.2 
we inclu de I Zw 18 as one of our sources (see also [Schruba et| 
[ar]|2012| ). The open triangles for |Schruba et al.| ( |2012| ) corre^ 
spond to their compilation of literature data and the filled ones 
to their metal-poor dwarf observations. CO non-detections are 
shown by upward arrows. The two horizontal dashed lines cor¬ 
respond to the regressions with unit slope shown in Fig. and 
the dotted one to a robust regression as described in the text: 
log(SFR/L^Q) = (-2.25 + 0.15) [12-rlog(0/H)] +(11.31 + 1.3). 


5. Molecular depletion times, the CO-to-H 2 
conversion factor, and specific SFR 

Although we have shown that is correlated with SFR even 
at low metallicity, the main parameter of interest is H 2 mass, 
Mh 2 - In fact, the study of low-metallicity star formation and 
H 2 content is hampered by the unknown factor, aco^ nsed to 
convert observed CO intensities to Mh 2 . Observationally, Mh 2 , 
Qfco, and the time required to consume the molecular gas reser¬ 
voir, Tdep, are intimately related and difficult to disentangle. 
By definition, Tdep depends inversely on the star-formation rate 
Tdep = Mh 2 /SFR = ofco/SFR, but the potential simplicity of 
predicting Mh 2 from SFR (or vice versa) through Tdep is com¬ 
promised by the necessary CO-to-H 2 factor. The observed (or 
directly derivable) quantities are and SFR; thus only the ra¬ 
tio, TdepMco. is constrained. 

Much effort has been devoted to finding the dependence 


of aco on metal abundance (e.g., Wilson 

19951 [Israel[[1997 

Roso 

lowsky et al.|2003 Bolatto et al.|2008 

Leroy et al.|2009b 

2011 

[Genzel et al.[[2012t [Schruba et al.[[2012l [Bolatto et al.[ 

2013 

1 . High-resolution cloud virial mass 

measurements tend 


[Rosolowsky et al.|2003[|Bolatto et al.|2008| ). However, because 


CO is found only in dense regions, rather than in diffuse gas, 
aco depends on the physical scales probed by th e observations, 
especially at low metallicity ( Rubio et~ST]|1993| ). Given the rel¬ 


atively con stant Tdep in n earby star-forming spirals (e.g., Bigiel 
|et al.|200^ |Leroy et ST]|2008[ |Bigiel et SI]|2Q11| ), some groups 
have assumed a constant Tdep to infer aco in metal-poor galax¬ 
ies. Such an assumption tends to give a steep dependence of aco 
on Q/ H: aco~^ ( [Genzel et al.j2012 Schruba et al. 


2012| ). An alternative approach relies on IR observations to de¬ 


rive Mh 2 and infer aco by estimating the dust-to-gas ratio and 
assuming it is the same in adjacent regions with/without H 2 . This 
method has the considerable advantage of being independent of 
CO observations altogether, but requires comparable quantities 
of both Hi and H 2 , and relies on the assumption that emissivity 
properties of dust grains are the same in dense and diffuse envi- 
ronme nts. Recent work using this method, introduced by [Israel 


( 1997| ), also gives super-linear slopes, although somewhat shal¬ 
lower (and with large scatter) than those obtained with_the_^ 
sumption of constant Tdep: (^co~^ oc(Z/ZQy ^-^~^ (Gratier et al. 
2010[|Leroy et al.|2011[ Bolatto et al.|2011| ). 


Here we present a different approach to constrain aco at low 
metallici ty. We exploit the results by |Saintonge et al. ( |2011b| ) 
(see also [Huang & Kauffmann||2014| ) who find a significant in¬ 
verse correlation of Tdep with sSFR: galaxies with higher sSFR 
also have shorter molecular gas depletion times. However, sSFR 


tends to als o be inversely correlated with met allicity (e.g., [Salim 
et al.|2Q14] ), as shown in Fig. [^(adapted from [Hunt et ^.|2Q1^ 


Thus, we must find a way to disentangle the effects of sSFR on 
Tdep from those of metallicity. 
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Fig. 6. Specific SFR vs. 12 +log(0/H) for our sample and includ¬ 
ing data from the literature ([Hunt et al.[2012t[Schruba et al.[2012[ 


[Filho et al.|2Q13[[BothweH et al.|2014 ). The metallicities for the 

ALLSMOG ( Bothwell et al.[2014[ ) saniple have been m odified as 
discussed in the text. As in Fig.[4[ the open triangles for Schruba 
et al. (2012 [ ) correspond to their compilation of data from the 


literature and the filled ones to their metal-poor dwarf observa¬ 
tions. The regression line is a robust best-fit of the Local Volume 


Legacy (LVL, [L^e et al.[2009 1 and the dwarf irregulars from 1 __ 

[et al.[ ( [2006[ ). The various colors for the LVL give the Hubble 
type (early types red, magenta; late types blue, cyan). 
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Hunt et al.: Molecular depletion times at low metallicity 


5.1. The CO-to-H 2 conversion factor 


[Saintonge et al.| ( |2011b ) found that star-forming galaxies of 
lower mass have much shorter Tdep than more massive ones, 
while the depletion time for atomic gas remains relatively con¬ 
stant. They also found that r^ep is strongly correlated with sSFR, 
with a single re lation being able t o fit both “normal” (“main- 
sequence”, e.g., Salim et ^|2007| ) galaxies and more extreme 
starbursts including LIRGs and ULIRGs. Their relation is based 
on a single value of aco^ appropriate for the properties of their 
sample, selected by stellar mass. 



Fig. 7. L^q/SFR vs. sSFR for our previous IRAM s ample and 
including data from the literature (CQLDGASS: [Huang & 
Kauffmann||2014| [Gao & Solomon||2004t |Garcia-Burillo et al. 


2012t|Schruba et al.|2012t|Hunt et al.|20i4 ). Our best-fit of the 

COLDGASS-rGao&Solomon-rGarcia-Burillo samples is shown 
(which is consistent with the regression given by |Saintonge et al. 
( 2011b| ) for only the COLDGASS galaxies). 


Here we explore the hypothesis that the Saintonge et al. cor¬ 
relation between Tdep and sSFR extends to the lower-mass (and 
higher sSFR) galaxies in our sample. However, for a better com- 
parison with our SFR estimat es, we use the COLDGASS data by 
[Huang & Kauffmann ( 2014| ) who re-evaluated the SFRs within 


the IRAM beam using additional mid-infrared and UV data for 
COLDGASS galaxies. Because there is some evidence that the 
SFRs from fitting of spectral-energy distributions such as done 



^ We have performed the entire analysis with both data sets and con¬ 
fir m an ove restimate of -0.2 dex in the SFRs meas ured by | Saintonge et| 


( |2011b| ) relative to [Huang & Kauffmann| < |2014] ). This change affects 
only slightly the metallicity depende nce of aco, steepening the p ower- 
law index by 0.08 with respect to the Huang & Kauffmann ( 2014| ) data. 


also consider other data from the literature where both CO mea¬ 
surements and stellar masses are available (see Sect. |2.5| ). The 
dashed line in Fig. [7] indicates our fit with a robust linear regres¬ 
sion to CO LDGASS data (using the SFR s and Mstar within a 22" 
apertu re by [Huang & Kauffmann|2014|), together with the sam - 


ples of Gao & Solomon (2004|) and Garcia-Burillo et al. (2012]): 
log(L^Q/SFR) = (-0.38 + 0.03) log(sSFR )-r( 4.75 + 0.31). The 
new fit is consistent with that found by [Huang & Kauffmann 
( [2014[ ) for only the COLDGASS galaxie s. However, both the 
dwarf galaxies in our sample and those in [Schruba et al. ( 2012 ) 
show large residuals relative to the general trend. 

The residuals to this fit are plotted vs. sSFR and 0/H in Fig. 
M Because Tdep = o^co L^q/SFR, the residuals from the corre- 
mtion of L^q/SFR on sSFR should eliminate the dependence 
on Tdep: [ L^q/SFR ](observed)/Tdep(fit) = [ Q'coCZVq'coCZq) ]"^ 
Thus, these residuals should correspond directly to the inverse 
of aco, relative to Z©. The left panel of Fig. [^demonstrates 
that there is no additional correlation with sSFR not accounted 
for by the COLDGASS relation; however, the right panel shows 
that there is a clear correlation with 0/H. The correlation of the 
residuals with metallicity is able to almost completely rectify 
the deviation of the low-metallicity galaxies by almost two or¬ 
ders of magnitude from the COLDGASS galaxies. We find a 
robust fitted slope of aco~^ ^ (Z/Z©)^-^^-^-^^, more shallow than 


found with the assumption of constant Tdep ([Schruba et al. 2012 


Genzel et al.|2Q12[ ), but consistent with the aco metallicity de 


pendence found by using dust continuum measurements ([Leroy 
[et al.|2011l[Bolatto et al.|20~TT] ). 

For the calculation of Mh 2 , we adopt a solar value 
aco = 3.2 M© (Kkm s“^ pc^)~\ not including a factor of 1.36 for 
helium (see also [Saintonge et al.[2011b[). For our galax i es and 
those from [Schruba et al.[ ( [2012 ) and [Bothwell et ^ ( [2014| ) 
having metallicities Z/Z©<1 (12-rlog(0/H)<8.69), we applied 
aco^(^/^Q)~^'^^', foi* metallicities XIZq>\, we used the con¬ 
stant aco as given above. Figure [^ shows the result of converting 
our L^q observations to Mh 2 - Although there is some scatter, the 


trend of Tdep with sSFR esta blished by [Saintonge et al.[([2011b|); 
Huang & Kauffmann ( [2014 ) seems to extend to the higher sSFRs 


and shorter Tdep of these dwarf galaxies. A steeper dependence 
of aco on metallicity would be inconsistent with this trend, but 
more data at 12-rlog(0/H)^8 are needed to better assess any sys- 
tematics. 

In any case, in agreement with previous work (e.g.. 


Saintonge et al. 

[201 Ibl [Bothwell et al. 20141 [Huang & 

Kauffmann 

2014 

), our analysis shows a strong variation of 

the molecular gas depletion times, with values ranging from 


-lOGyr for the COLDGASS galaxies to ~70Myr for IIZw40 
and even smaller, <10 Myr for SBS 033 5-052E, NGC 5253 and 
NGC 1569 (from [Schruba et aL][2012[ ). This variation is corre¬ 
lated with a range in sSFRs of almost three orders of magnitude. 

The right panel of Fig. [^suggests that Tdep depends weakly 
on metallicity, even after me aco dependence has been ac¬ 
counted for. The regression line is the robust fit to our galaxies 
and those from [Schruba et al.] ( [2012| ); the correlation is signifi¬ 
cant (Pearson coefficient of 0.38, significance >99.9%), with a 
slope of 0.64+0.11. A similar (weak) effect was also found by 
Leroy et ( [2013a[ ), but over a narrower range in metallicities 
and Tdep. Because Tdep and sSFR are correlated, and sSFR and 
0/H are also correlated (see Fig.[^, such a result is perhaps not 
surprising; it could be the result of the secondary dependence on 
sSFR or on Mstar- In Sect. [^ we discuss forms of aco and its 
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Fig. 8. Residuals from the regression shown in Fig. 


vs. sSFR (left panel) and vs. 12+log(0/H) (right). We include here data for 


our sample and literature data (CQ LDGASS: [Huang & Kauffmann|20T4]|Gao & Solomon|2004[[Garcia-Burillo et al.|2()T2 ; |SchruM 


|et al. 12012} |Bothwell et al.||2014| ). The horizontal dashed line in the left panel shows 0 to aid the eye, while the regression in the 


right panel is the robust fit to the significant trend of the residuals of the L' 


CO- 


/SFR fit with metallicity cx(Z/Zq)^-^^ 



sSFR [yr-*] 12+log(0/H) 


Fig. 9. Tdep vs. sSFR (left panel) and vs. 12+log(Q/H) (right) for our sample and including data from the literature (CQLD GASS: 


|Huang & Kauffmann||2014t |Gao & Solomon||2004t 

Garcia-Burillo et al.||2012t |Schruba et al.||2012 Bothwell et al.||20141). As in 

previous figures, the open triangles for|Schruba et al. 

(|2012|) correspond to their compilation of data from the literature and the filled 


ones to their metal-poor dwarf observations. The regression line is our fit to the |Schruba et al.|(2012| ) and our sample as described 
in the text. 


metallicity dependence (e.g., |Wolfire et al.||2010t |Krumholz et| 
|al.|2011D . 

Because 0/H and Mstar are correlated through the mass- 
metallicity relation, in Fig. 1^ we have investigated the possibil¬ 
ity that Tdep depends on Mstar- This correlation is also formally 
significant, even more so than the one with 0/H (Pearson coef¬ 
ficient of 0.56, significance >99.999%). The slope of 0.46+0.03 
(shown by a dashed line in Fig.[^ is similar to that (0.41+0.05) 
found for the COLDGASS galaxies alone. Our results suggest 
that low-mass metal-poor galaxies tend to have high sSFR (Fig. 


and shorter H 2 depletion times (Figs.[^[^, and the two prop¬ 
erties are correlated. This result is not altogether new, although 
we have shown that it extends to much lower stellar masses and 
higher sSFRs than the COLDGASS sample. 

6. Atomic and molecular gas content 

The implication is that for low-mass, metal-poor galaxies, with 
high sSFR, Tdep is systematically shorter than for more mas¬ 
sive, metal-rich systems. Because our sample combined with 
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Fig. 10. Tdep VS. Mstar for otir sample and including data from 
the literature ([Huang & Kaulfmann||2014t [Schruba et al.||2Q12} 


Fig. 12. Mgas = (Mhi + MH2)/Mbar vonic VS. Mgtar, including data 


Bothwell et al.|2014|). A s in previous figures, the open triangles 

for |Schruba et al.| ( |2Q12| ) correspond to their compilation of data 
from the literature and the filled ones to their metal-poor dwarf 
observations. Our robust best-fit regression for these samples is 
shown as a dashed line; the dotted line corresponds to the slope 
for the COLDGASS galaxies alone ( [Huang & Kaulfmann|2014[ ). 


from the literature (COLDGASS: [Huang & Kauffmann 


the [Schruba et al.[p012[ ) and the Filho et aL] ( 2013[ ) compilation 
spans a wide range of 0/H and sSFR, with ^^CO(l-O) detections 
down to Z/Zq~ 0.1, it is well suited to assess the dependence on 
these parameters of molecular and total gas (Hi-rH 2 ) content, 
relative to stellar mass and total baryonic mass. The H i ma sses, 
Mhi, for otir sample are giv en in Table [T] (s ee also Sect. [2.4[ ) and 
we have take n Mhi for the [Schruba et ah] ( [2012[ ) galaxies from 
Walter etH^ (TH INGS, [2008| ) when available, or otherwise from 
Kennicutt et al.[ ( [2Q0^ . 

We first investigate any trends of MH 2 -to-total gas mass frac¬ 
tion in Fig. both as a function of sSFR (left panel) and 
12-hlog(0/H) (right). We use total Hi mass for this comparison, 
although it is almost certainly more extended than the H 2 (see 
Table [T]), which usually resides within the optical confines of the 
galaxy ( [Regan et~ar|[2001[ [Leroy et ST|[2009a[ ). The combined 
samples suggest a weak trend for H 2 fractions to decrease at low 
sSFR, but there is no trend with 0/H. A better assessment of 
both trends would require more CO data from high-sSFR and 
metal-poor galaxies (and resolved H i measurements). 

Fig. [T2| shows the total gas Mgas (= Mh 2 +Mhi) fraction rela¬ 
tive to the total baryonic mass Mbaryonic (= Mgas+Mstar) as a func¬ 
tion of stellar mass. As in Fig. m we use the total H i mass for 
our targets (Col. 12 in Table [T]). By including also the data from 
Schruba et al.[(2012]);[Filho et al.[ ( [2013] ); [Bothwell et al.[ ( |2014[ ); 
Huang & Kauffmann| ( 2014[ ), we find there is a significant cor¬ 
relation between gas mass fraction and Mstar, as shown by the 
robust regression line. This is mainly due to the spread in Mstar 
of the two most extreme samples, COLDGASS and [Filho et al.[ 
( [2013[ ). At low gas-mass fractions, this may also be partially in¬ 
duced by the co-dependence of the two axes, because of the pres¬ 
ence of Mstar- Nevertheless, the trend is roughly continuous over 
a factor of 10^ of Mstar, from ~ 10^ to ~ 10^^ Mq. 


[Schruba et al. [20121 [Filho ot al. [201 3 [[Bothwell et al.[2Q14| ). Our 

robust best-fit regression for these samples is shown as a dashed 
line. 


Two of our sample galaxies are clear outliers, Mrk 996 and 
NGC7077, both of which have a smooth elliptical-type mor¬ 
phology despite their low metallicity. Thus they may have al¬ 
ready exhausted much of their H i, since the H 2 -to-total gas mass 
fractions are relatively high, -40%. These two galaxies could be 
in the process of quenching their star formation, unless more gas 
is accreted from the outside. 

We examine in Fig. the trend of total gas mass rela¬ 
tive to Mbaryonic VS. sSFR (left panel) and 12-i-log(0/H) (right). 
Here also the total gas fraction is significantly correlated with 
both sSFR and 12-rlog(0/H) (see the robust regression lines), al¬ 
though there are significant outliers in our sample for the metal¬ 
licity trend, which is mainly dominated by the XMP galaxies. 
However, as in Fig.[^ there may be a co-dependence of the two 
axes because of the presence of Mstar (left panel) and because of 
the stellar mass-metallicity correlation (right). 

The Hi depletion times, Tdep(Hi), for our combined sample 
are shown in Fig. There is more than an order-of-magnitude 
scatter, but no systematic trend relative to sSFR. This result 
is similar to that found for more massive galaxies, and the 
horizontal dashed line in Fig. corresponds to the roughly 
constant Tdep (Hi) o f 3.4 G yr found for the GASS galaxies by 
Schiminovich et al. ( 2010[ ). Fig. [Td] shows a possible weak trend 


with 0/H, which starts to affect Tdep(H i) at 12-rlog(0/H):^8. This 
could be a consequence of the lar ge fraction of H i-dominated 
galaxies at these metallicities (e.g., [Krumholz|20f3) ). 

In general, we find that galaxies are more gas rich at low stel¬ 
lar masses, at high sSFRs, and at low metallicities, although the 
metallicity trends ten d to have more scatter. Similar conclusions 
have b een reached by Bothwell et al. ( [2014| ) and by [Filho et al. 
( 2013[ ), although for a more limited range in parameter space. 
Our combined samples span five orders of magnitude in Mstar, a 
factor of - 10^ in sSFR, and a factor of -100 in oxygen abun¬ 
dance. The trends for total gas fraction are similar to those we 
find for the molecular gas depletion time Tdep: low-mass, metal- 
poor, and high sSFR galaxies tend to have shorter Tdep- Tdep de- 
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sSFR [yr-i] 12+log(0/H) 


Fig. 11. MH 2 /Mgas = Mh 2 /(Mhi + Mh 2 ) vs. sSFR (left panel) and vs. 12+log(0/H) (right) for our sampl e and including data from th e 
literature ( [Schruba et al.|2012t|Bothwell et al.|2014[); the left pane l gives also the COLDGASS data ( Huang & Kaulfmann|2014| ). 
As in previous figures, the open triangles for |Schruba et al.| ( |2012} correspond to their compilation of data from the literature and 
the filled ones to their metal-poor dwarf observations. 



sSFR [yr-i] 12+log(0/H) 


Fig. 13. Mgas = (Mhi + MH2)/Mbary 

onic VS. sSFR (left panel) and vs. 12-rlog(0/H) (right) for our sample and including data from 

the literature (Schruba et al. 

2012 

|Filho et al.||2013t|Bothwell et al. 

2014]); the left panel gives also the COLDGASS detections 

(Saintonge et al.|2011bl). As in previous figures, the open triangles for 

Schruba et al. (|2012|) correspond to their compilation of data 


from the literature and the filled ones to their metal-poor dwarf observations. The regression lines shown in both panels correspond 
to our fit of the galaxies from Schruba et al.| ( |2012 ); Filho et ar] ( [2Q13| ), ALLSMOG, and our sample as described in the text (with 
all the uncertainties alluded to in Sect. |2.5| ). 


creases by a factor of -100 over a variation of -1000 in sSFR 
and a factor of -10 in 0/H. However, Hi depletion times are 
roughly constant, although with large scatter, across this entire 
range of sSFR and 0/H. 


7. Discussion 

The short molecular gas depletion times, Tdep^lOOMyr, we find 
at the highest sSFRs (^0.1 Gyr“^) are comparable to starbursts 
at z ^ 1; they are smaller th an those found for “ main-sequence” 
galaxies at similar redshifts ( [Genzel et al.|2015| ). Thus, the low- 
metallicity galaxies studied here may be low-mass analogues of 
starbursts observed at cosmological distances. This is not un¬ 
likely since most of them are likely merging or interacting sys- 
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sSFR [yr-‘] 12+log(0/H) 


Fig.l4.Tdep(Hl) vs. sSFR (left panel) and vs. 12+log(0/H) (right) for our sa mple and including da ta from the literature ( [Schruba et| 
al.|2012t Bothwell et al.|2014| ). As in previous figures, the open triangles for Schruba et al. ( 2012| ) correspond to their compilation 


of data from the literature and the filled ones to their metal-poor dwarf observations. The constant Tdep(Hi), shown by a horizontal 
dashed line, corresponds to 3.4 Gyr as found for the GASS sample ( [Schiminovich et al.|2010| ). 


terns (see Sect. |2.1| ). A starburst could enhance star-formation 
efficiency, namely the quantity of H 2 converted into stars over 
a free-fall time (e.g., [Krumholz et '^^12009 1; mergers or inter¬ 
actions would be expected to drive gas into the central regions, 
causing higher gas densities which could, at least partly, com¬ 
pensate for the low metallicity through increased self shielding. 

Independently of the origin of the short molecular gas de¬ 
pletion times, the secondary dependence of Tdep on Mstar and/or 
metallicity is not straightforward to interpret. Because of the way 
we have derived the molecular gas depletion times, the determi¬ 
nation of a power-law dependence of aco on metallicity is in¬ 
timately related to Tdep. To better understand whether short Tdep 
and its dependence on Mstar and 0/H are artefacts of our deriva¬ 
tion of aco, here we investigate an alternative formulation of 
aco, and its potential relati on with Tdep- Wh at we have defined 
as Tdep is the inverse of what | Krumholz et al.|(|201 If call specific 
gas SFRs, SSFRgas(CO,H 2 ,total), the ratio of SFR surface den¬ 
sities to the various gas components of the ISM: Egas, and 
Ehi. Given a total gas mass surface density, Egas, Krumholz et al. 
give precise predictions for the fraction of molecular gas mass 
in the form of H 2 ,/h 2 , as well as the fraction of molecular gas 
where CO dominates over other forms of carbon,/co- 

In their mod els (see also [Krumholz et al.|[200^ [McKee &| 
Krumholz|2010|), /h 2 is a function of the ratio of the interstellar 


radiation field (ISRF) at the surface of the gas cloud, Gq, and the 
hydrogen number volume density, Hh, as well as the dust optical 
depth to UV photons, Tdust- /h 2 depends on metallicity through 
Tdust which is assumed to be linearly proportional to Z/Zq. The 
transition in a cloud of atomic and ionized carbon to CO is driven 
mainly by dust extinction, thus cloud optical depth. The fraction 
of gas mass where the carbon is in CO form,/co, depends on/H 2 , 
modified by an exponential function which depends only weakly 
on Gq, and hh, but strongly on metallicity through Ay (i.e., Tdust). 


Th us, by a ssuming a characteristic value for G'Jhh, Krumholz et 
aLjpOlip derive the dependence of the gas SSFRs on gas mass 
surface density, Egas, and metallicity, Z/Zq. 

Fig. shows their predictions together with the data 
we have analyzed here. Although the models are formu¬ 
lated in terms of surface densities, to better compare with 

( |2011[ ) include 


global values 
correction factor 


Krumholz et al 


( 0 , 


a resolution 

which we have adopted according to 


their prescription. Following [Krumholz et al.[ ( [2011| ), the data 
points correspond to molecular gas masses, Mh 2 (CO), in¬ 
ferred from CO luminosities using a constant conversion fac¬ 
tor, aco = 4.4 M© (Kkm s“^ pc^)“^; this is the same aco as used 
up to now (see [Saintonge et al.|2011b| ) but here with the inclu¬ 
sion of helium. The CO-inferred molecular mass has been de¬ 
rived assuming a total gas surface density Egas = 20 M© pc“^; as 
shown by the blue curves in F ig.p?) the exact value o f Egas is 
not critical. As emphasized by Krumholz et al.[ ( [2011 ), this in¬ 
ferred mass Mh 2 (CO) should not be interpreted as the total mass 
of molecular gas, but rather the mass of the molecular gas where 
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the pre dominant form of carb on is CO. Also shown in Fig 
are the Krumholz et al. (2011) model prediction^ as a function 
of metallicity, Z/Z©, for three values of tota l gas column density 
(for more details see [Krumholz et al.|2011| ). The models of CO- 
inferred H 2 mass differ from the true H 2 mass by a multiplicative 
factor that varies exponentially with Z/Z© (because of the form 
of/co described above, see also[Wolfire et aL][2010t[Bolatto et| 
al.[[2013[). As in [Krumholz et al.[ ([2011[), the SSFRgas(H 2 , CO) 
curves are shown only for metallicities and total gas surface den¬ 
sities that are sufficiently high that the H 2 and CO mass fractions 
are non-zero. Like the data compilation presented by [Krumholz 
eTaTld^OTT] ), our data are also qualitatively consistent with the 
trends predicted by the models. 


Wolfire[ 


assumed to vary linearly with Z/Z© (for more details see[ _ 

et al. 1201 Ob. In these models, SFR surface density depends on/H 2 
and Sgas, and is defined by a roughly constant timescale, 2.6 Gyr. et al. ( 2011| ). 


^ He re we continue to use the 8.69 solar calibration of [Asplund et al. 
([2009 1 rathe r than the value of 12+log(0/H) = 8.79 used by [Krumholz 
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Fig. 16. Ratio of SFR and CQ-inferred H 2 mass vs. relative metallicity as described in the text, i ncluding our data and data from 
the literatur e (CQLDGASS : Huang & Kaulfmann|2014t[Schruba et al.|2012t[Bothwell et al.|2014[ ). As in previous figures, the open 
triangles for Schruba et al.| ( |2012| ) correspond to their compilation of data from the literature and the filled ones to their metal-poor 
dwarf observations. In the left panel, Mh 2 is inferred with our formulation of ago as a power - law dependence on Q/H as des cribed 
in the text; in the right, we use the exponentially-d ependent ago as found by Wolfire et al.| ( |201()| ); |Krumholz et al. ( 2011| ). Also 
shown are the predictions by '' . . 


Krumholz et al. 


(2011), showing three values for total gas surface density, Egas: 5 M© pc 20 M© pc 


and SOM© pc indicated by dashed, dotted, and solid lines, respectively. Red color shows total gas, purple H 2 , and blue CO. 



Z/Zo 


Figure shows that the data and the models agree ex¬ 
tremely well for the prediction of CO luminosity; however at 
significantly sub-solar metallicities, Z/Z©^0.2, there are some 
discrepancies relative to the H 2 prediction, with some data ex¬ 
ceeding by a factor of 3 or more the expected variation of the ago 
factor. This is shown in Fig. 16 where we compare the power- 


law dependence of ago derived previously with the exponen¬ 
tial variation predicted by the models. Unlike in Fig, p?] where 
Mh 2 (CO) depended on a constant ago, here we derive Mh 2 from 
a metallicity-dependent ago ; the left panel shows the power-law 
dependency ago ocZIZq~^-^ as shown in previous fi gures, and the 
right p anel the exponential depend ency found by Wolfire et al. 
( 2010| ) and |Krumholz et al.] ( |2Qll| ), agooc exp(-Z/Z©). As dis¬ 
cussed above, such a dependence arises through the assumption 
that ago depends exponentially on the visual extinction within 
the clouds. Ay, and that Ay is proportional to Z/Z©. 


Fig. 15. Ratio of SFR and CO-inferred H 2 mass vs. relative 
metallicity as described in the text, including our data and data 
from the literature (COLDGASS:|Schruba et al.|2012l|Bothwell| 
^TaLlIMT^. As in previous figures, the open triangles for 
Schruba et al.| ( |2Q12| ) correspond to their compilation of data 
from the literature and the filled ones to their metal-poor dwarf 
observations. Following [Krumholz et al.| ( |2011| ), here Mh 2 is 
inferred with a constant ago =4.4M© (Kkms“^ pc^)“\ as de¬ 
scribed in the text. Also shown are the predictions by |Krumholz| 
et al.| ( |2011| ), showing three values for total gas surface den¬ 
sity, Sgas! 5M©pc“^, 20M©pc“^, and 80M©pc“^ indicated by 
dashed, dotted, and solid lines, respectively. Total gas corre¬ 
sponds to red, purple to H 2 , and blue to CO. 


Because both Mh 2 and SFR predicted by the Krumholz et al 


(2011) models depend onfn^, the expectation is that SSFR(H 2 ) 
should be constant, independently of metallicity. Fig, shows 
that the data do not quantitatively follow this trend, but rather 
show a residual increase of SSFR(H 2 ) toward low metallicities. 
Neither formulation (power-law or exponential) of ago gives a 
horizontal trend of SSFR(H 2 ) as a function of Z/Z©. The expo¬ 
nential ago seems to better approximate the models around so¬ 
lar metallicity, but overestimates ago at Z/Z©~0.1. On the other 
hand, the power-law ago reproduces well the lowest metallic¬ 
ity SSFR(H 2 ) (perhaps by construction, given that our ago was 
heavily infiuenced by the lowest metallicity galaxy in our sam¬ 
ple). In any case, it seems that the short depletion times, corre¬ 
sponding to high SSFR(gas) are not a result of our approach, nor 
is our power-law formulation of ago significantly inferior to the 
exponential one. 
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As shown in Fig.[^ galaxies with XIZq^Q.I are not always 
consistent with the model predictions for SSFR(H 2 ). Because 
of the large span of Tdep [inverse of SSFR(H 2 )] in our sample, 
it could be that the roughly constant timescale in the models 
is an oversimplification; indeed the Krumholz models are not 
designed for starbursts with potentially short depletion times 
and a high-density ISM. Moreover, there is strong evidence that 
clumping factors ar e different for molecular and atomic phases 
( Leroy et al.|2013b| ), implying that the use of a single correction 
factor to compare the predictions by |Krumholz et al.| ( |2Ql 1| ) with 
observations is almost certainly an oversimplification. 

There could also be other problems with the applicability of 
the model predictions to these metal-poor starbursts, since aco 
could have additional parameter dependencies beyond metallic¬ 
ity. The assumpt ion of constant G'Jhh m ay not be valid at low 


metallicities; the Krumholz et al. (2011) models assume Gq = 1, 
but the average U V tlux Gq at low metallicities is expected to be 
harder and, at a given total stellar mass, more intense because of 
the warmer stellar temperatures. There could also be some prob¬ 
lem with the assumption that dust opacities and the capacity of 
the gas to self-shield from UV radiation vary linearly with metal¬ 
licity. As mentioned in Sect. even at a similar low metallicity, 
galaxies can have dramatically diff e rent dust-to-gas and du st-to- 
stars mass ratios ( [Hunt et al.|2014| ). [Schneider et al.| ( |2015| ) sug¬ 
gest that this difference is due to the differences in the density of 
the cool gas in the ISM, because of the increased efficiency of 
grain formation in high-density environments. Thus, in addition 
to metallicity, ISM volume density could also be an important 
parameter for governing both the Hi-to-H 2 and the H 2 -to-CO 
transitions in galaxies. However, the incidence of high ISM den¬ 
sities at low metallicities is not yet well established. More obser¬ 
vations of the molecular ISM are needed for metal-poor galaxies, 
especially at Z/Zq^ 0.1 -0.2 where the validity of CO as a tracer 
of H 2 begins to be called into question. 


8. Summary and conclusions 

We have detected ^^CO(l-O) with IRAM 30-m observations of 
a sample of 8 metal-poor dwarf galaxies, including a tentative 
detection at a metallicity of 12-hlog(0/H) = 7.7, ~0.1 Z/Zq. By 
calculating stellar masses and SFRs for our sample, and com¬ 
paring our observations with a compilation of data from the lit¬ 
erature, we have ex tended the trend for mole cular Tdep with stel¬ 
lar specific SFR by Saintonge et al.| ( |2011b| ) to lower metallici¬ 
ties, lower stellar masses, and higher sSFRs. In agreement with 
[Saintonge et al.| ( [2()IIb| ) and [Huang & Kauffmannj ( [2QT4) , we 
conclude that Tdep is not constant in galaxies; in our combined 
sample Tdep varies by a factor of 200 or more (from ^50 Myr to 
~10 Gyr) over a range of 10^ in stellar sSFR and Mstar- 

By exploiting the correlation of Tdep with sSFR, we were able 
to constrain aco as a function of metallicity, and found a depen¬ 
dence of aco ^ similar to that found by work ba sed on 

dust-continuum measurements compared with gas mass ( [Leroy[ 
et al.|20ir|[Bolatto et al.|20lT] ). In conjunction with atomic gas 


measurements taken from the literature, we used our formulation 
for aco to infer total gas masses, and total baryonic masses. The 
gas mass relative to the total baryonic mass depends on Mstar, 
stellar sSFR, and 0/H; we found significant correlations for all 
three parameters in the sense that gas-mass fractions are higher 
for galaxies of lower stellar mass, lower metallicity and higher 
sSFR. However, despite large scatter, the H i SFR depletion time 


is roughly constant for our combined sample, ~3.4 Gyr (see also 
[Schiminovich et al.|2010[ ). 

We have exami ned our data in the c ontext of the star- 
formation models by Krumholz et al. ( [2011[ ), and compared the 
power-law formulation of aco with the theoretically motivated 
exponential formulation. The predicted ratios of SFR and gas 
mass compare extremely well for the CO data, but the SSFR(H 2 ) 
inferred from either form of aco are inconsistent for Z/Zq-^0.2. 
A residual trend of SSFR(H 2 ) remains even when it should be 
invariant with metallicity according to the models. The cause of 
the disagreement is not clear but could arise from the starburst 
nature of our metal-poor sample; more observations of molecu¬ 
lar tracers at low metallicity are needed to perform a more ex¬ 
tended comparison with model predictions. 
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